home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / ctbcon.z / ctbcon
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCTTTTBBBBCCCCOOOONNNN((((3333FFFF))))                                                          CCCCTTTTBBBBCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CTBCON - estimate the reciprocal of the condition number of a triangular
  10.      band matrix A, in either the 1-norm or the infinity-norm
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CTBCON( NORM, UPLO, DIAG, N, KD, AB, LDAB, RCOND, WORK, RWORK,
  14.                         INFO )
  15.  
  16.          CHARACTER      DIAG, NORM, UPLO
  17.  
  18.          INTEGER        INFO, KD, LDAB, N
  19.  
  20.          REAL           RCOND
  21.  
  22.          REAL           RWORK( * )
  23.  
  24.          COMPLEX        AB( LDAB, * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      CTBCON estimates the reciprocal of the condition number of a triangular
  28.      band matrix A, in either the 1-norm or the infinity-norm.
  29.  
  30.      The norm of A is computed and an estimate is obtained for norm(inv(A)),
  31.      then the reciprocal of the condition number is computed as
  32.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  33.  
  34.  
  35. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  36.      NORM    (input) CHARACTER*1
  37.              Specifies whether the 1-norm condition number or the infinity-
  38.              norm condition number is required:
  39.              = '1' or 'O':  1-norm;
  40.              = 'I':         Infinity-norm.
  41.  
  42.      UPLO    (input) CHARACTER*1
  43.              = 'U':  A is upper triangular;
  44.              = 'L':  A is lower triangular.
  45.  
  46.      DIAG    (input) CHARACTER*1
  47.              = 'N':  A is non-unit triangular;
  48.              = 'U':  A is unit triangular.
  49.  
  50.      N       (input) INTEGER
  51.              The order of the matrix A.  N >= 0.
  52.  
  53.      KD      (input) INTEGER
  54.              The number of superdiagonals or subdiagonals of the triangular
  55.              band matrix A.  KD >= 0.
  56.  
  57.      AB      (input) COMPLEX array, dimension (LDAB,N)
  58.              The upper or lower triangular band matrix A, stored in the first
  59.              kd+1 rows of the array. The j-th column of A is stored in the j-
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCTTTTBBBBCCCCOOOONNNN((((3333FFFF))))                                                          CCCCTTTTBBBBCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              th column of the array AB as follows:  if UPLO = 'U', AB(kd+1+i-
  75.              j,j) = A(i,j) for max(1,j-kd)<=i<=j; if UPLO = 'L', AB(1+i-j,j)
  76.              = A(i,j) for j<=i<=min(n,j+kd).  If DIAG = 'U', the diagonal
  77.              elements of A are not referenced and are assumed to be 1.
  78.  
  79.      LDAB    (input) INTEGER
  80.              The leading dimension of the array AB.  LDAB >= KD+1.
  81.  
  82.      RCOND   (output) REAL
  83.              The reciprocal of the condition number of the matrix A, computed
  84.              as RCOND = 1/(norm(A) * norm(inv(A))).
  85.  
  86.      WORK    (workspace) COMPLEX array, dimension (2*N)
  87.  
  88.      RWORK   (workspace) REAL array, dimension (N)
  89.  
  90.      INFO    (output) INTEGER
  91.              = 0:  successful exit
  92.              < 0:  if INFO = -i, the i-th argument had an illegal value
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.